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We explore the connections between the constraints on the precision of quantum logical operations 
that arise from a conservation law, and those arising from quantum field fluctuations. We show that 
the conservation-law based constraints apply in a number of situations of experimental interest, such 
as Raman excitations, and atoms in free space interacting with the multimode vacuum. We also 
show that for these systems, and for states with a sufficiently large photon number, the conservation- 
law based constraint represents an ultimate limit closely related to the fluctuations in the quantum 
field phase. 



I. INTRODUCTION. 

It was pointed out in 0] that, when trying to do quan- 
tum logic with a quantized field, three possible kinds of 
errors could arise: errors due to qubit-field entanglement 
(a point also made in 2]), to field phase fluctuations, and 
to field amplitude fluctuations. The last two are easy to 
separate from each other, but it is not immediately clear 
how they relate to the first. 

Almost simultaneously, and quite independently, it was 
also pointed out 3] that the existence of a conservation 
law, which holds for many of these quantized-field sys- 
tems, would also constraint the accuracy achievable in 
some quantum logical gates. It is natural to ask the ques- 
tion of how this result may relate to the ones derived in 
by very different methods, and this is, broadly speak- 
ing, the purpose of the present paper. 

The remainder of this introductory section is devoted 
to some preliminary observations, followed by a brief 
summary of the rest of the paper's sections and results. 

It seems that the conservation-law induced quantum 
limit (which we shall abbreviate as CQL from now on, 
for consistency with is most directly related to the 
entanglement error: if the conservation law is of the form 



a) a = const 



(1) 



then it is impossible, starting from a factorized state 
that is an eigenstate of a z , to produce a pure eigen- 
state of <j x — that is, to accomplish a transformation like 
l°> (1°) + l 1 ))/^ without getting entangled with the 
field. Intuitively, it is clear that if the transition from |0) 
to |1) requires the absorption of a photon, then the state 
of the field associated with the state 1 1) in the above su- 
perposition will, in general, be different from the state 
associated with |0). Nonetheless, the actual amount of 



entanglement depends on the field state and the coupling 
assumed, and the question has proved to be somewhat 
subtle (in particular, for coherent field states, which are 
not changed by the loss of one photon). We refer the 
readers to the work of Silberfarb and Deutsch Q on this 
question, and also to the papers 0, H] for other related 
aspects of the controversy. (We will comment later on 
another issue that was clarified in 0,13], namely, the re- 
lationship between the constraints in |l| and those posed 
by spontaneous emission for an atom in free space.) 

There is still the question of how this relates to phase 
and/or amplitude fluctuations. The form of the CQL is 
(in a way that we shall make precise later; see Section II) 
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4 1 + Ao-{nf 



(2) 
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where o~{n) is the standard deviation of the number op- 
erator n = a^a in the initial state of the field. P e is the 
"error probability" for the quantum logical operation, de- 
fined as the largest deviation (in the norm-square sense) 
between the intended final state of the qubit and the ac- 
tual final state, maximized over all input qubit states. 
Eq. J5J indicates that this error probability decreases as 
the fluctuations in the field intensity increase. Thus, it 
seems clear that this error is not due to field-intensity 
(or, probably, field-amplitude) fluctuations. There is, 
however, a suggestion that (at least in some limit) it 
might be related to phase fluctuations, since for these 
a relationship of the form a((f>) 2 > l/4a(n) 2 is gener- 
ally understood to hold for some reasonable approximate 
Hermitian phase operator <f> (see, for instance, and 
also the discussion below, in Section IV). 

The suggestion is, therefore, that the CQL J5J may be 
related to the error caused by fluctuations in the field 
quantity conjugate to the one that appears in the con- 
servation law l|T[l. This is explored, and a precise cor- 
respondence is established, in Section IV of this paper. 
Prior to this, in Section II, we present a simple proof of 
the result and in Section III we show some examples 
of systems of interest where a conservation-law-induced 
constraint applies. Of particular interest are a Raman- 
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coupled system (where the conserved quantity is less ob- 
viously related to the total angular momentum), and an 
atom in free space, where spontaneous emission is the 
principal limiting factor. Section V briefly summarizes 
our conclusions. 



II. DERIVATION OF THE CQL. 

This section presents a simplified derivation of the con- 
straint on quantum logic arising from a conservation law, 
specialized for our present purposes; see ^(| f° r a rnuch 
more detailed derivation and additional context. 

In what follows we assume that we are trying to per- 
form a quantum logical operation (typically a Hadamard 
gate) on a qubit using a Hamiltonian that involves the 
qubit as well as a quantized control field, subject to a 
conservation law of the form L = L\ + L 2 = const, with 
Li a qubit operator and L 2 a field operator. Typically 
L\ will be the Pauli operator a z , written in the compu- 
tational basis {|0), 1 1)} as 



<7 Z = |0){0|-|1)(1| 



(3) 



The identification of the states |0) and |1) with appro- 
priate states of the qubit (typically an atom) will depend 
on the system (see next Section for examples). 

We therefore assume that the field-qubit interaction 
Hamiltonian is such that after a certain time t the evo- 
lution operator U (t) (the time argument will be omit- 
ted below for simplicity, as long as no ambiguity results) 
approximately carries out a Hadamard transformation, 
while leaving invariant the quantity 



L = a z + L 2 



(4) 



We note that the exact Hadamard transform goes as fol- 
lows: 



1 



|0>--^(|0> + |1» 
|1>-^(|0>-|1» 



(5) 



and therefore it changes the operator a z into a x — 
|0)(1| + |1)(0|. Accordingly, we may take as a measure of 
the success of the gate the smallness of the difference (or 
"error" ) operator D defined as 



D = U*a z U - <j x 



(6) 



Consider now the commutator \<j x + D,L] = [U'cr z U, L]. 
We have 

[a x + D,L} = U f a z UL - LU ] a z U 

= U^a z UU^LU - U^LUU^a z U 
= tf[o- z ,L]U 

= (7) 



since, by assumption, LU — L, and a z commutes with 
L (it clearly commutes with itself, and it also commutes 
with L 2 , at equal times, because they act on different 
systems). From this we conclude 



[L,D] = [a x ,L}=-2ia y 



(8) 



and therefore the standard deviations of D and L in an 
arbitrary (initial) state of the joint qubit-field system 
must satisfy 

a(D)a(L)>l\{[a x ,L})\ = \{a y )\ (9) 



The right-hand side has its maximum value (equal to 1) 
in an eigenstate of a y , so we find 



a(Df T 



> 



1 



1 



a(L) 2 l + a(L 2 ) 



(10) 



where the last equality follows if one assumes that the 
initial state of field and qubit is a product state (note 
a(cr z ) 2 — 1 in an eigenstate of a y ). 

Equation (|f 0[l is the basic constraint on the accuracy of 
a Hadamard gate in the presence of a conserved quantity 
of the form @. In the remainder of the paper we shall 
work mostly with this form, but it is useful to relate 
the error measure cr(D) 2 nax to the more familiar "gate 
fidelity" F that is defined in terms of the overlap of the 
final qubit state with the intended state. As 1 — F 2 gives 
the probability to find the qubit in a state orthogonal 
to the intended one, this quantity may be termed the 
failure (or "error") probability of the gate, P e . Taking F 
to be the "worst case" fidelity, over the set of all possible 
initial states of the qubit, it can be shown (see Appendix 
A below) , along the lines of [lj} > that 



P e = 1 — F 2 > -(D 2 



1 



1 2 1 

As an example, for a system for which the quantity Q 
(proportional to the total angular momentum of the sys- 
tem) is conserved, we would write L = a z ± 2a'a > so 
L 2 = 2aJa, and then (jl 1|) would immediately yield the 
result ©■ We discuss this further, and provide several 
additional examples, in the following Section. 

Although the discussion above has been very specific, 
it is straightforward to see how all these arguments can 
be generalized to arbitrary systems and operations. The 
key ingredients are: (1) a conserved quantity L that can 
be written as the sum of a qubit operator L\ and a con- 
trol operator L 2 , and (2) a desired evolution that would 
change L\ into another qubit operator that does not com- 
mute with Li(0). One then finds that (3) such an op- 
eration cannot, in general, be carried out exactly, and 
the unavoidable error (for at least some initial states) 
is inversely proportional to the sum of the uncertainty 
squared of L\ and L 2 in the initial state. 
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III. EXAMPLES OF SYSTEMS WHERE THE 
CQL APPLIES. 

A. Two-level atom and single-mode field 
(Jaynes-Cummings model). 

The simplest system where a conserved quantity of 
the form 0) is found to hold is the well-known Jaynes- 
Cummings model ^lj' which describes the coupling of a 
two-level atom with a single-mode of the electromagnetic 
field. Allowing for a detuning A between the atom and 
the field, the Hamiltonian for the system may be written 
in a suitable interaction picture as 

H = hAafa + ihg (a® \e){g\ - a) ® \g){e\) (12) 

where g is an appropriate coupling constant. To see how 
this Hamiltonian performs an approximate Hadamard 
gate, replace the field operators a and a* by classical 
constants Eqc 1 ^ and i?oe _I ^, respectively, and make the 
identification |e) — > |1), \g) — > |0). Then the evolution 
operator, for zero detuning, is 



U{t) = cos(gE t) 







-i<p 



o 



sm(gE t) (13) 



in the {|0), |1)} basis, and at the time t 
4> = 0, this becomes 



n/4gE , if 



(14) 



This is not actually a Hadamard gate H, but rather a x H 
(the state |1) is changed to |1) — |0)). Nonetheless, it is 
easily verified that the operation U^a z U produces — er x , 
and hence the results in the previous section apply if the 
difference operator D is defined as D = U'a z U + a x and 
the sign of a x is changed in the equations that follow . 

It is trivial to verify that the quantity —\<J Z + a) a 
is conserved under the Hamiltonian Ijl2(l (regardless of 
whether the detuning A is zero or not), and therefore 
Eq. (QUI holds, with L 2 = -20*0. By Eq. JTTJ, this 
yields a failure probability 



P P > 



4 1 + 4cr(n) 2 



(15) 



which means that when the field is quantized it is im- 
possible to implement the transformation (|14|) exactly. 
Note the similarity between Eq. i|15|) and the contribu- 
tion of phase fluctuations to the error probability esti- 
mated in Eq. (8) of Q, for a state with the optimum 
<r{(t>) 2 = l/4cr(n) 2 . This is elaborated further in Section 
IV below. 



B. Two-level atom and multimode field. 

After the publication of [lj , Itano [fj called into ques- 
tion the applicability of the results obtained from an "ef- 
fective single mode" model to the important situation of 



an atom in free space, and argued that for such a system 
spontaneous emission would be the only source of error 
for quantum logic, provided the field was in a coherent 
state. Although in a certain sense Itano's claim regard- 
ing spontaneous emission is accurate, it was shown in Q 
that the constraints obtained from the Jaynes-Cummings 
model still apply. This can be seen easily in the present 
formalism by considering the Hamiltonian 



H = H 



k 



Akaj^k + ih ^2 9k ( ak <S> |e) (g\ -a J k (g> \g) (e 

k 



(16) 

for an atom interacting with many quantized field modes, 
labeled by the index k, through coupling constants <7k- 
This can describe, in particular, an atom in free space 
interacting with a laser field, and allowed to emit spon- 
taneously into the available vacuum modes. 

It is straightforward to verify that for the Hamiltonian 
(|16fl . with the same identification |e) — > |1), \g) — ► |0) as 
before, the following quantity is conserved 



-L 



2 4«k 



(17) 



and therefore a constraint of the form (|llfl holds, with 
L-2 = — 2J2k n k- If the initial field state is a product 
state over all the modes, then oiL-i) 2 — 4 J^w ^("k) 2 - If, 
moreover, all the modes are in a coherent state (of which 
the vacuum is a special case), then cr(nk) 2 = n.k, where 
flk is the average number of photons in mode k, so the 
constraint takes the form 



1 — F > 



1 1 



4 1 + An 



(18) 



where n is the total number of photons in the field. This 
can be directly compared with other multimode estimates 
such as [lj| and pjj l. 

In his comment |fj , Itano pointed out a seemingly un- 
physical aspect of the result i|18|) , namely, the fact that 
one might apparently improve the accuracy of the opera- 
tion just by making the laser beam wider, while keeping 
its intensity constant: this would formally increase the 
number of photons in the field, and hence decrease the 
right-hand-side of l|18|). but in fact it would not change 
at all the evolution of the atomic system. This point is 
certainly true, but it does not mean that Eq. (|18|) is an 
incorrect lower bound, only that it is not (in general) a 
very tight lower bound. As shown in Q, when the effects 
of spontaneous emission are properly taken into consid- 
eration, as Itano suggested, one finds a constraint on the 
fidelity that can be put in the same form as (|18fl . only 
with a generally smaller number of photons n', of the 
order of the number of photons actually intercepted by 
the atom's scattering cross section. That is to say, one 
typically has 



1-F J > 



1 



> 



1 1 



4 1 + 4n 



(19) 
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While the question of how to achieve the lower bound 
(|18fl will not be taken up in detail until the following Sec- 
tion, we may already point out that there is, clearly, a 
lot of information in the Hamiltonian (|16[) that does not 
show up anywhere in the conservation law or in the con- 
straint (|18|l : this includes, in particular, all the coupling 
constants t/k, which account for the experimental geome- 
try and, ultimately, for whether the coupling of the atom 
to the field is optimal or not. It is natural to expect that 
(|18|l will only be a tight lower bound on 1 — F 2 when the 
coupling between the field and the atom is as close to op- 
timal as possible. In free space, for an atom interacting 
with a paraxial beam, the coupling is far from optimal 
and there are a lot of "wasted" photons; hence it is not 
surprising that a tighter lower limit on 1 — F 2 should be, 
under those circumstances, rather larger than the one ob- 
tained from the conservation law alone. Note, however, 
that one could, in principle, account for this suboptimal 
coupling in a somewhat "ad hoc" way by modifying the 
conserved quantity (|T7ll so as to include only the pho- 
tons in those modes that effectively couple to the atom, 
by introducing an appropriate "effective mode decom- 
position" that is better adapted to the geometry of the 
atomic radiation pattern (for instance, one based on vec- 
tor spherical harmonics, as opposed to plane waves). 



C. Raman-coupled three-level atom. 

A popular arrangement for quantum logic involves a 
three-level atom in which leves \g±) and |e) are coupled 
by a dipole transition (and, experimentally, by a laser 
field of suitable frequency and polarization), and levels 
|<ji2) and |e) are also coupled by a dipole transition (and 
a laser field of their own). The qubit states would then 
be \gi) = |0) and \g 2 ) = |1), so that now 



IffiXffil - I02X52I 



(20) 



The simplest Hamiltonian for this system involves two 
field modes, a and b, each with a possible detuning, 

H =hAi a' a 

+ ihg a (a <g> |e) (g 1 1 - a f <g> \gi) (e|) 

+ ihg b (b®\e)(g 2 \-tf ®\g 2 )(e\) (21) 

One has now, for the commutator [a z , H] , 

[a z ,H] = - ihg a (a <g> |e) (gi\ + a f <g> |si)<e|) 

+ ihg b (b(g> \e)(g 2 \+tf ®\g 2 )(e\) (22) 

whereas 

[a)a,H] = -ihg a (a® |e)(fln| -I- a f ® \gi)(e\) 
[tfb,B\ = -iUg b (6® \e){g 2 \ +tf ® | fl2 )(e|) (23) 

and hence it is clear that the appropriate conserved quan- 
tity is 



(Note that this could have been obtained from conser- 
vation of angular momentum considerations, if, for in- 
stance, the states \gi) and \g 2 ) corresponded to magnetic 
quantum numbers m = 1 and m = — 1; then the pho- 
tons in modes a and b must have opposite circular po- 
larizations, and hence opposite angular momentum com- 
ponents along the quantization axis.) Using again 
with L — b^b ~ a^a and the assumption that the modes 
are initially independent, one obtains the constraint 



1 - F 2 > - 



4 l + <j(n a ) 2 + er(n b ) 2 



or, for coherent states 



1 — F > 



1 



1 



4 1 + n a + n b 



(25) 



(26) 



This generalizes immediately to a multimode, free- 
space case, along the same lines as in the previous sub- 
section, by introducing modes at that couple to the 
I <7i ) — * |e) transition, and modes 6k that couple to the 
\g 2 ) — * |e) transition; the fact that \g\) and \g 2 ) are or- 
thogonal implies that one can choose the modes so as to 
form two disjoint sets (e.g., they would have orthogonal 
polarizations). The result would again be of the form 
()25p. only with a sum over modes in the denominator of 
the right-hand side. 



IV. SATURATING THE CQL. 

As pointed out in the previous Section, the CQL may 
be well below the actual achievable error probability for 
many systems. This is because many systems (i.e., many 
different Hamiltonians) may, in fact, satisfy the same con- 
servation law, and the CQL derived from it must hold for 
the worst as well as for the best coupling possible. The 
purpose of this Section is to determine what this "best" 
possible coupling is, in the sense of the one that gets 
as close as possible, in practice, to the error probability 
illl . In doing so, we shall naturally learn what are the 
basic physical properties of the system that are captured 
by the CQL (fTT). 

To begin with, as pointed out in the previous sec- 
tion, we shall want to have optimal coupling, meaning a 
single-mode (or "single-effective-mode" ) situation; hence 
we start by working with the Jaynes-Cummings Hamil- 
tonian H12|l and asking how close it actually gets to the 
limit H15|) . The question can be addressed by noticing 
that the full evolution operator for the Hamiltonian l|12|) 
can be written, again through direct exponentiation of 
H, as 



L = a. 



tfb- 



(24) 



U 




(27) 
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from which follows 

D = U^<T z U + a x 

I cos 2gtVaJa 
\ 1 — sin 2gtV 'aa+ —i= 



1 — sin 2gt\> at a 




— cos 2gW aa 



Now, we want the dispersion cr(D) 2 = (D 2 ) — (D) 2 of D 
in a state such as |+y)|^>), where \+y) = (|0)+i|l))/\/2 is 
an eigenstate of a y (for which the right-hand side of 10 is 
maximized), and \ip) is an arbitrary field state. Actually, 
we shall write 



|0>=^C„|n) 



(29) 



with real coefficients C n , since we got the ideal transfor- 
mation H14|) by letting the classical field be real. In this 
case we find that 

1 



(D) = - ((cos2gt^/n) - (cos 2gWn + 1)) (30) 

For a state such as a coherent state with a large average 
number of photons n, we expect that (|3Ull will be of the 
order of 1/n, since it involves the difference between y/n 
and yjn + 1, which is of the order of n -1 / 2 , multiplied 
by gt, and we expect gt ~ w/i^/E at the time when the 
approximate transformation l|14(l is realized. Then, when 
we square H30|) to compute a(D) 2 , we will obtain some- 
thing of the order of 1/n 2 , which for large n is negligible 
versus the leading terms in l|15(l (of the order of 1/n). 
Hence, we will neglect this term in what follows. 

For the square of D, we obtain, from Eq. l(2*%]l. the 
simple result 

D 2 =2 — s'm2gt\ // a^a =a) — sin 2gt\/ aot =a 

va^a vaot 

+ a x ^cos 2gtVaUi — cos 2g£v / aaT^ (31) 
and thus we have, in the state | + y)\i/j}, 



(D 2 ) = 2-2^2 C n G n+ i sin(2.gi\/^TT) (32) 



Naturally, by the normalization of \ip), we must have 
^C 2 = 1. Now we may consider the various reasons 
why the second term on the right-hand side of 132|) might 
not be exactly equal to 1 around the time t — n/Agy/E. 
One reason is the spread in Rabi frequencies, g\fn\ for 
different values of n, it is impossible to make all the sines 
in the sum (KS2[) equal to 1 at the same time. Another 
reason is that we are adding the products C„C„+i, rather 
than C 2 . The first reason clearly means that the system 
is sensitive to intensity fluctuations (uncertainty in the 
photon number). The second reason, as we shall show 
below, captures the system's sensitivity to phase fluctu- 
ations. 

To estimate how much of a(D) 2 , in the Jaynes- 
Cummings model, is due to each of these two sources of 



error, we can write C n C n+ i = C 2 + {C n C n+ \ — C 2 ), and 
sin(2gt Vn+T ) = sin^tV" + 1) + {sm(2gt^/n + 1) - 
sm(2gty/fi + 1)), and expand (|32[l. neglecting the prod- 
uct of terms expected to be small. The result is 

<j{D) 2 ~2 - 2 s\n(2gty/n+l) ^ C n C n+ i 

n 

-2^2 Cl(sm(2gWn + l) - sm{2gWn + 1)) 

n 

(33) 



Now let 2gty/n + 1 = 7r/2. 
approximated by writing n 



The second sum can be 
n + An, and expanding 



sm(2gt\Jn + 1) in a power series in An/(n + 1), to sec- 
ond order; the sum ^2C 2 (An) 2 then equals tr(n) 2 . The 
result is 



a(D) 2 ~2-2^C n C n+1 



16n 2 



(34) 



The first sum, on the other hand, can be expressed in 
terms of the Susskind-Glogower "exponential of phase" 
operators 0], exp(i(j)) = \n)(n + 1|, as (exp(i</>)). As- 
suming that the phase of the field is a reasonably well- 
defined quantity with zero average, this can be expanded 
as 



^C„C„ + i = (exp(i(/))) ~ 1 



= 1 



a{4>) 2 



We thus get the final result 



a{Df 



16n : 



: a(n) 2 + a{ct>f 



(35) 



(36) 



for the contribution of the field intensity and phase fluc- 
tuations to the gate error in the Jaynes-Cummings model 
(compare this to the result of a semiclassical calculation 
in Eq. (8) of 0). Assuming the approximate validity of 
c(0) 2 > l/4<r(n) 2 (which, for a coherent state with large 
n, holds as an approximate equality), we get 



a{Df 



> 



Mnf 



4cr(n) 



(37) 



and it is clear that the second term in l|37|l is very 
much like the conservation-law constraint I|1UI) . with 
L-2 = (±)2n, as long as one can assume cr(n) 2 ^> 1. 
Again, this is automatically the case for a large-n co- 
herent state, and, at any rate, it is the most interesting 
case, since one typically will want the QCL to be very 
small, and this can only happen if cr(n) 2 3> 1. 

On the other hand, there is clearly nothing like the first 
term in (|36|l in the conservation-law constraint l|l(J|) ; that 
is, there is no indication of an error that grows with the 
size of the intensity fluctuations. This is simply because 
the atom-field coupling that results from the Hamiltonian 
111 21 has an intrinsic photon-number dependence that is 
not a necessary consequence of the conservation law . 
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One can, in fact, envision a whole family of Hamiltonians 
compatible with the same conservation law, of the general 
form 

H = ihgJ2(fn\n)(n + 1| ® \e)(g\ - f n \n + l)(n\ ® | 5 )(e|) 

n 

(38) 

where the /„ are arbitrary numbers. The Jaynes- 
Cummings model corresponds to the choice /„ = \/n + 1, 
but one could equally well (formally, at least) remove 
all dependence on the photon number by setting /„ = 
1 for all n. Recalling again the Susskind-Glogower 
cxponential-of-phase operator, we see that the resulting 
Hamiltonian is just what one would obtain by quantiz- 
ing only the field phase in the semiclassical Hamiltonian 
H = ikgE (e' l ' li \e}(g\ — e~' l ^\g}(e\) discussed in connection 
with Eq. (JT3J): 



H = ihg^{\n){n 



II 



\e){g\-\n+l){n\®\g){e\) 



(39) 

We could continue the discussion with the Hamiltonian 
(|39[) (see Appendix B for details), but it is more transpar- 
ent at this point to work with a hypothetical Hermitian 
phase operator that satisfies 

[kn] = ~i (40) 

and rewrite our (also hypothetical) interaction Hamilto- 
nian (|39ll as 



H 



ihg(e^\e){g\ 



-i<p 



\9)(e\ 



(41) 



(for clarity, in the remainder of this section the operators 
4> and n have been given a carat). 

Apart from making the correspondence with the semi- 
classical case more obvious, the form (|4*T|) has the sat- 
isfactory feature of making the conservation law (JTJ 

straightforward, since by l|4(J|) the operators e ±l< ^ are "dis- 
placement operators" for the photon number h, changing 
it by one unit up or down. Thus, in a sense, Eq. (|41|l is 
the simplest nontrivial Hamiltonian compatible with a 
conservation law of the form QJ: it just captures the es- 
sential fact that the photon number must change by ±1 
as the atom makes a transition up or down. 

It is known (again, see 9] for details) that there are 
difficulties with the postulated Hermitian operator (f>, es- 
sentially, because one cannot lower the number of pho- 
tons below the vacuum level "|0)." Pegg and Barnett 
have defined a Hermitian <fr for which something like l|4(J|) 
approximately holds for most relevant, "physical" states 
[l5| . One of us ^(| has also shown that it is possible 
to define a Hermitian phase operator in a larger Hilbcrt 
space, augmented by "macroscopic" states with a (for- 
mally) infinite photon number, using nonstandard anal- 
ysis. We shall not go into these complications here, but 
merely assume that something like Eqs. (|4(JH and l|41|) 
makes sense, in an appropriate limit, for useful, physi- 
cal states (as the derivation in Appendix B shows, it is 



enough to assume that the contribution of the vacuum 
state to the field is negligible). 

The evolution operator associated to the interaction 
Hamiltonian (|41|l is then just like the semiclassical result 
(|13[) . only with a quantized <f>: 



U(t) = cos(gt) 




sm(gt) (42) 



and thus at the time t = n/Ag it becomes 



U 



__L K - 

V2 [e^ 



—i<f> 



(43) 



In an initial state | -I- y)\ijj), the difference operator D = 
U'a z U + <J X has the expectation value (D) = (sin0), 
which we can set equal to zero, and then we find 



a(D) 2 = (D 2 ) = 2 - 2(cos< 



(44) 



To minimize o~(D) 2 , we clearly want the phase distribu- 
tion for the field state to be centered at zero and 
as narrow as possible. Under this assumption of small 
phase dispersion (which clearly implies large o~(h) 2 , by 
100} ), we can expand l|4~4l and get 



a(D) 2 



> 



2—2(1— 

a{^) 2 
1 




4a (n) 2 



(45) 



We have, therefore, finally succeeded at devising a sit- 
uation where the lower limit <r(D) 2 > 1/(1 + 4a(n) 2 ) can 
be approximately reached (assuming large a(n) 2 ). The 
necessary ingredients were: (1) an effective single-mode 
interaction; (2) an atom-field coupling that is insensitive 
to intensity (i.e., photon number) fluctuations, and (3) 
a state with minimum uncertainty in the phase-intensity 
variables. In the large n (and small <j(4>) 2 ) limit, the 
CQL clearly captures the error due to (quantum) field 
phase noise. 



V. CONCLUSIONS 

In this paper, we have presented a number of systems, 
involving the interaction of an atom with the quantized 
electromagnetic field, for which constraints on the accu- 
racy of quantum logical operations can be derived from 
conservation laws of the form first discussed by one of us 
in |3j. These constraints would appear to include those 
arising from atom-field entanglement (although we have 
not explored this isssue in detail), and, in the appropri- 
ate limit, they seem to include as well those arising from 
fluctuations in the field phase. The special significance of 
the phase is clear from the "minimal Hamiltonian" l|41|) 
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compatible with the conservation law as the conju- 
gate variable to the photon number, which is the field 
quantity that appears in the conservation law, it must 
be involved in the Hamiltonian in the way given by (|4H , 
to ensure that photons are created and destroyed appro- 
priately as the atom makes a transition. This leads to 
its unavoidable presence in the evolution operator, as in 
Eq. i|43[l , which then naturally becomes sensitive to phase 
fluctuations. 

It follows from the above that quantum logical oper- 
ations that are not constrained by a CQL of the form 
considered here must be (at least to lowest order) insen- 
sitive to fluctuations in the field phase. 

On the other hand, we cannot conclude that operations 
not constrained by the CQL can be performed without 
any errors, since the possibility still remains that they 
might be sensitive to intensity fluctuations, whose effects, 
as shown in the previous section, are not accounted for 
by the CQL. Whether this is the case or not, and if so, to 
what extent, clearly needs to be ascertained separately 
for each specific physical system under consideration; the 
present work merely provides a starting point by showing 
why such an undertaking would be necessary. (See, in 
this context, the discussion in 0,0.) 

Finally, it is not yet, at this point, entirely clear how 
these considerations may relate to the ultimate limits to 
"quantum logic with quantum controls" discussed by one 
of us in • As elaborated in later work, it seemed that 
an important ingredient in the derivation of such lim- 
its was the requirement that the interaction should be 
"turned on and off" by the quantized control system it- 
self. It would be interesting to see how much of that work 
might be clarified by the present formalism; for instance, 
by formulating some of the constraints in terms of other 
conservation laws, by searching for approximate conser- 
vation laws, or by separating the constraints that follow 
from such laws from those that have a different origin. 
We expect to be able to present results on this program 
in the near future. 

The authors acknowledge support from the U. S. Army 
Research Office, the Strategic Information and Commu- 
nications R&D Promotion Scheme of the MPHPT of 
Japan, the CREST project of the JST, and the Grant- 
in- Aid for Scientific Research of the JSPS. 



where the field states \Eg) are not necessarily normalized. 
Next, assuming that the initial qubit state is | + y) (as 
we did in the derivation of Eq. (|1 C)|) ) we have 

(D 2 ) = ^\(+y\(U^ z U-a x ) 2 \+y)\^) 

= \\(U^ z U-a x )\+yM\\ 2 

= \\(a z U-Ua x )\+y)\m 2 

= ||| J B °)-| J B 1 )|| 2 + |||i??) + | J B 1 1 )|| 2 (A2) 

where the last equality follows from the direct application 
of (|A1|) to the states | +y)\ip) and a x \ +y)\^). 

Now, using again IjAlfl . observe that the first term of 
(|A2(I can be written as 

|||£; )-|i? 1 )|| 2 = 2||(0|[/|-x)|V)|| 2 (A3) 

where | — x) = (|0) — \l))/^/2 is the —1 eigenstate of 
a x ; hence this term expresses the probability that the 
qubit may be found in state |0) at the end of the trans- 
formation, if the initial state was | — x). Note, however, 
that a pure Hadamard gate should turn | — x) into |1) 
exactly; hence this term is nothing other than twice the 
error probability (or "infidelity"), 1 — F 2 (| — x)), for an 
initial state | — x) . Similarly it can be seen that the other 
term in l|A2[) equals 2(1 — F 2 (| + x))). We conclude that 

{D 2 ) ^ A~2F 2 {\~x))~2F 2 {\+x)) <A{l~F 2 mn ) (A4) 

when the minimum fidelity is taken over all the possible 
initial states of the qubit. (Note that there is a typo- 
graphical error in the corresponding equations (88) and 
(90) in 10]: both fidelity terms must subtract from 1, as 
in Eq. (|A4I) above.) 

Finally, noting that (D 2 ) > cr(D) 2 always holds, we 
obtain 

1-FL»> \{D 2 ) >\o-{D) 2 (A5) 
which is Eq. (fTTf 1 - 



APPENDIX B: DERIVATION OF (JUJ) WITHOUT 
ASSUMING A HERMITIAN PHASE OPERATOR 



APPENDIX A: DERIVATION OF THE FIDELITY 
EQUATION (fTTT) 

The purpose of this Appendix is to show in detail how 
one can go from the noise cr(D) 2 of Eq. I|10|) to the "error 
probability" 1 — F 2 of Eq. ijTT|) . The derivation closely 
follows 10], circa Eq. (88). First, the action of the evolu- 
tion operator U in the computational basis can generally 
be written as 

U\0)\tl>) = \0)\E%) + |1>|E°) 

C/|1)|V) = |0)|^) + |1)|^) (Al) 



The evolution Hamiltonian H39|l can be easily inte- 
grated to yield the evolution operator 

U = cos gt 

( 2sin 2 (. 9 i/2)|0)<0| -smgtj: n \n + l)(n\\ 
\singtY, n \n){n + l\ J 

(Bl) 

In the expression above, the state |0) represents the pho- 
ton vacuum. The effect of this extra term can be made 
negligible if one works with a field state that has a negligi- 
ble projection on the vacuum; for instance, for a coherent 
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state, Co — e n l 2 . Neglecting it altogether, then, we see 
that av is transformed as 



U^aM = 



cos2gt — sixx2gt^2 n \n + l)(n 

■ sin 2gt^ n \ n){n + 1| —cos2gt 

(B2) 



Without any further approximations, one then obtains, 
for the expectation value of I? = U^azU +<t x , in the state 

(D) =-ism(2 g t)(^2(\n + l)(n\-\n)(n+l\)^ (B3) 

which equals zero if the state \tp) = ^ n C n \n) is chosen 
to have all real coefficients C n (zero average phase) . For 



the operator D 2 , one gets 

D 2 =2 - sm(2gt) ^(|n + l)(n\ + \n)(n + 1|) 

'1 0' 



sin'(2 5 i)|0)(0| <g> 







(B4) 



Thus, neglecting again the field's vacuum-state contribu- 
tion, we get for cr{D) 2 



a{D) 2 = 2-2^2 C n C n+ i sin(2<7*) = a{4>f 



> 



1 



4a(n) 2 
(B5) 

at the time gt = 7r/4, where the connection with phase 
fluctuations is established via the Susskind-Glogower op- 
erators, as in the derivation of l|35|) . 
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